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Abstract. In this paper we analyse both the dynamics and the high density physics of the infinite di- 
mensional lattice gas model for random heteropolymers recently introduced in fll . Restricting ourselves to 
site-disordered heteropolymers, we derive exact closed deterministic evolution equations for a suitable set 
of dynamic order parameters (in the thermodynamic limit), and use these to study the dynamics of the sys- 
tem for different choices of the monomer polarity parameters. We also study the equilibrium properties of 
the system in the high density limit, which leads to a phase diagram exhibiting transitions between swollen 
states, compact states, and regions with partial compactification. Our results find excellent verification in 
numerical simulations, and have a natural and appealing interpretation in terms of real heteropolymers. 

PACS. 61.41.-|-e Polymers, elastomers and plastics - 75.10.Nr Spin-glass and other random models 



1 Introduction 

The practical advantage of infinite dimensional models 
for random heteropolymers over the more standard ones 
(which are either defined on finite dimensional lattices, or 
in real space) is the absence of the polymer chain con- 
straint. In infinite dimensional (or 'lattice gas') models, 
where every site is defined to be a neighbouring site of 
every other site, the constituent monomers of a polymer 
are given unrestricted freedom to occupy positions. The 
movements of the latter are controlled only by energetic 
considerations, by the maximum number of monomers al- 
lowed to occupy any given site, and by trivial restrictions 
on the number of moves allowed to occur simultaneously. 
In a recent paper jl], which was devoted to the equilib- 
rium solution of a particular infinite dimensional model, 
this property was shown to enable an exact solution, in 
spite of the presence of quenched disorder, (viz. the real- 
ization of the constituent monomers of the chain, which 
are assumed to come in different species. The quenched 
disorder is entirely specified by the proportions of the var- 
ious monomeric species present. No further aspects of the 
system are quenched). In those cases where one's interest 
is mostly in issues which are not directly related to spa- 
tial properties of the heteropolymer, the compromise of 
infinite dimensionality might be a price worth paying for 
exact solvability. 

The present paper is concerned with the further anal- 
ysis and exploration of the model introduced in to 
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which we refer for a more detailed and extensive discus- 
sion of the relevant literature and of the motivation un- 
derlying the study of random heteropolymers in general 
(e.g. their role as precursors of models for proteins), and 
lattice gas models of random heteropolymers in particu- 
lar (see also e.g. ^,§,D and references therein). We would 
like to emphasize that this model does not represent re- 
alistically the polymeric features of a heteropolymer, but 
it does capture some of the characteristics, and for this 
reason it can be considered as a precursor for a model 
of random heteropolymers. We complement the equilib- 
rium study in two distinct ways. Firstly, we solve the 
dynamics of the model (which we define as a stochastic 
process of the Glauber type), by deriving exact closed de- 
terministic laws in the thermodynamics limit for a suit- 
ably chosen set of dynamic order parameters. This can 
be done, at least for site-disordered models, as a direct 
consequence of the simplification induced by the infinite 
dimensionality of the model, and allows for an exact quan- 
titative analysis of the approach to equilibrium. Secondly, 
we study the equilibrium properties of the model of p| 
in the limit of high local density. Since the original proce- 
dure followed in Q for solving the model (which applies to 
both site-disordered and bond-disordered systems), does 
not allow for this limit to be taken, we first present an 
alternative (but equivalent) equilibrium solution (which 
applies only to site disorder) where taking the high local 
density limit is allowed. We then find, after an appropriate 
re-scaling of the system temperature (in order to prevent 
entropy dominance of the free energy) a non-trivial phase 
diagram with a distinct swollen phase, a compact phase, 
and a phase characterised by partial compactification of 
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the heteropolymer. Especially near the transition lines in 
this diagram, the system is found to exhibit a highly non- 
trivial and slow plateau- type dynamics. 

We supplement our theoretical analysis by extensive 
numerical simulations, which support our predictions re- 
garding the dynamics and regarding the high-density phase 
diagram very satisfactorily (within the experimental lim- 
itations imposed by the slow dynamics). 



2 Microscopic Rules of the Model 

The model of describes a system of N monomers, la- 
beled i £ {1, --jN}, and R = aN possible sites for these 
monomers, labeled r S {1, ..,R}. The degrees of freedom 
of the system are the locations G {1, . . . , i?} of the N 
monomers, and the overall system state is defined by the 
state vector r = (ri, . . . , tat). Each site, however, can oc- 
cupy at most Uc monomers; this imposes R constraints on 
the state vector of the form J2i^ri,r < "-c (one for each 
r), and limits the global monomer density according to 
a > n~^. Those positions at a given site which are not 
occupied by monomers, are assumed to be occupied by 
solvent molecules. The monomers are allowed to interact 
(in pairs) only when they find themselves at the same site. 
Within equilibrium statistical mechanics this is described 
by a Hamiltonian of the form 



H{r) 



(1) 



The disorder in the system is in the realisation of the 
interaction energies Xij. Unlike j^, we will in this paper 
restrict ourselves to models with site disorder, where 



(2) 



such as the Random Hydrophobic Model (RHM) 
where 0(A^A') = ^{X + X'), or the Random Charge Model 
(RCM) |,||, where ^(A, A') = eX^Xj. Here the variable A; 
defines the monomer species of monomer z, with Ai > for 
hydrophylic/polar monomers, and Ai < for hydrophobic 
monomers. The RHM describes a system where hydropho- 
bic monomers try to avoid contact with the solvent, and 
all hydrophilic monomers try to optimise contact with the 
solvent. The parameter e e { — 1,1} in the RCM controls 
whether monomers of identical polarity attract or repel 
each other. We restrict ourselves, for simplicity, to the 
case whether the disorder is discrete in nature and where 
there are K monomer species, i.e. Xi G {Ai, ...,Ak}- 

We next endow this equilibrium system with a Glauber- 
type dynamics. Time is discretised, and the probability to 
find system state r at iteration step £ is written as Pi{r). 
The elementary transitions we will allow for are those 
where a single (randomly drawn) monomer is moved to 
a (randomly drawn) new site; to describe such transitions 
we define an operator Fi^ which moves monomer i to site 
x: Fix{r) = (ri, ...,ri_i,a;,ri+i, ...,r7v). The evolution of 
the system is a Markovian stochastic process of the form 



Pi+iir) = V W[r] r']pi{r') 



r 



(3) 



with the transition probabilities 

N R 



^n,(r'),«e \p{Fixr';r')5ri.x + [^-p{Fixr']r')\5^^y, | 



and 



p[r\r'] = 



e 2 



2cosh[i/3Z\i/(r;r' 



(5) 



5ij, and 



in which AH{r;r') = H{r) - H{r'), 5.,j = 1 
nx{f) = X^i ^ri.x (the number of monomers found at site 
x). The process (H) is one where we forbid all transitions 
which would lead to violation of the constraint n^, (r ) < ric, 
and where we allow all other single particle moves (where 
both the candidate particle i to be moved and the can- 
didate site X for it to move to are drawn uniformly at 
random) with Glauber- type probabilities, such that the in- 
variant distribution will be the Boltzmann state (provided 
the system is ergodic, which in the present case is intu- 
itively clear). One easily convinces oneself that the transi- 
tion matrix (^ is properly normalised, i.e. W[r, r'] = 
1 for all r', and that the Markov chain (^) indeed obeys 
detailed balance, i.e. W[r,r']pao{i'') — W[r' ,r]poo{r) for 
all {r,r'}, with Pooif) ~ exp[— /3iJ(r)]. 

Although the Markov chain (||) is the representation 
which is most conveniently implemented in numerical sim- 
ulations, for a theoretical analysis of the dynamics a con- 
tinuous time description would be preferable. This is achieved 
(sec JpO| ) upon choosing the duration of each iteration step 
of (^Pto be a Poisson distributed random number (with 
average r), which converts the Markov chain into a master 
equation of the form 



r^pdr) 



r' 



W[r;r']pt{r') - pt{r) 



Upon choosing r — N^^ (so that every monomer is tar- 
geted for update on average once per unit time), and upon 
inserting (0), we then arrive at the master equation 



ix r' 



(6) 



Using the properties of the Poisson distribution, one can 
show that the relative uncertainty on the time axis, in- 
duced by having random durations of iteration steps, will 
vanish for N ^ oo. 

We finally note that the only remaining non-trivial 
Glauber term (|^) occurring in (^ is of the following form: 

p[F„r; r]^^-^ tanh[i/3Z\„i/(r)] (7) 
with AixH{r) — H (Fi^r) — H (r) , which is here given by 
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Upon defining nx,xif) — X)i ^ri,x^\i,\ (counting the num- 
ber of type-A monomers at site x), we can write this as 

Z\„i/(r) - ^ 0(A„ A) 

X 

(8) 



3 Dynamics: Derivation of Closed 
Macroscopic Laws 

The appropriate dynamic order parameters of the system 
are the smallest set of macroscopic observables which obey 
closed and deterministic evolution equations in the ther- 
modynamic limit. In the present model this set turns out 
to be, at least for finite ric and for pair energies of the sep- 
arable form (H), the fraction cn{r) of sites with a given 
local population (or 'local state') n — {ha^, ■ ■ ■ j^yi/f ) of 
monomers. Here n\ indicates the number of monomers of 
species A. Thus 



cn{r) = 



1 



X—1 



Us 



.X=Ai 



(9) 



with n e {0, 



} . Note that J2n cn(^) = 1 and that 



cn(^) = as soon as J2\ "-a > "-c (due to our constraint 
on the maximum number of monomers at any given site) , 
for any r. We will write the constraint as g[n] S {0,1}, 
where g[n] = 1 if and only if J^x 't-a < "-c- 

The probability of finding the set of obervables {cnii")} 
at time t is given by 



(10) 



Its evolution in time follows upon taking a temporal deriva- 
tive and inserting the master equation m): 



r 



E 



l-tanh[i/3Z\„i?(r-) 



[ n 



Y\_S[cn-cn{r)] 
n 



We expand this expression, using 
CniFrxV) = cn(r) 



R 



X^Xi 



We also insert the partitionings 1 = Sn"[nA'^nA ^(r),n'^ 
and 1 = En"'[nA'5nA.rjr),n;;'], and we use (|). This then, 
together with the identity nx{r) = X^a "a,2;('"), leads to 



the following differential equation for Pt[{cn}]: 



d 
di 



men}] = -J:J- {Pdicn}] E E E -n"cn'" 
w " ^ V X n" n'" 

l-tanh[i/3[^ 0(A, A')«, -<;) + 0(A, A)]] 

A' 

X'^X . 
A'5^A J 

Taking the limit N —^ oo for finite ric and finite K con- 
verts this into a Liouville equation, describing determinis- 
tic flow, i.e. solutions of the form Pt[{cn}] = Hn'^I'''^ ~ 
cn(i)] where the deterministic trajectories {cn(t)} are 
simply given by the solution of the associated flow equa- 
tions. With the short-hands F^n = (nyi^ , . . . , n\-i, n\ ± 
1, TT-A+i, . . . , tiAk) these (coupled) flow equations can be 
simplified to 



A w 



rf^^" = EE^"' |%+nK + l)^^^„^,„/A[n'-n] 

+ c^-n"A^^^„,,„.+i£A 
- cnnx6y „, ex[F+n'-n] 



cnn'xS 



e^n-F^n'] (11) 



where 



^x[m] = i-itanh[i/3^0(A,A')mA'] 



(12) 



The imposed limit on the maximum number of monomers 
at a given site is seen to be built into (^l]) , since for any lo- 
cal monomer arrangement n with X)a ""A = "-c+l equation 
( ITTl ) indeed dictates that JjCn = if cn = 0. 

Since the dimensionality of our macroscopic laws in- 
creases with Tic and the number K of monomer types, fur- 
ther reductions can be helpful. One reduction is to con- 
sider the fraction p[n+,n_] of sites with local monomer 
arrangements where n+ monomers have Ai > and n_ 
monomers have Ai < (for the RHM this implies count- 
ing hydrophobic and hydrophilic monomers; for the ROM 
one counts monomers with a specific charge sign): 



p[n+, n_] - ^ cn (5y Sy ^ 



(13) 



with n± g {0,rtc}. We will, more specifically, inspect 



the (co-)variances 



(n±)^ and p = {n+n. 



They give information about the spatial properties of the 
system state: a large a± implies spatial concentration of 
a monomer type in a relatively small number of sites, 
whereas p informs us of the degree of spatial separation of 
the monomer types. These quantities will nor- 

mally not obey closed equations, but are calculated via 
( |l3|) from the (numerical) solution of (ITTl). 
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Fig. 1. Evolution of the order parameters {cn} for the Ran- 
dom Hydrophilic Model (RHM), with K — 2 (two monomer 
types, A € {—1, 1}, one hydrophobic and one hydrophilic), 
ric — 2, a = 1, and /3 = 100000. Upper graphs: theory, i.e. 
numerical solution of (|l^). Lower graphs: numerical simula- 
tions, with A*' — 3000. Note: in this example C(2^o) = for ''■U 
t. 

4 Dynamics: Comparison with Numerical 
Simulations 

We first compare the numerical solution of the macro- 
scopic laws ( IhI) with the results of measuring the observ- 
ables directly in numerical simulations of the micro- 
scopic process All simulations are presented as single 
runs, as oppose to ensemble averages. A typical example is 
shown in figure for the RHM (i.e. <j){X,X') = i(A + A'), 
with K = 2 (two monomer types, A G { — 1,1}, one hy- 
drophobic and one hydrophilic), Uc = 2, a = 1 and (3 = 
100000. In this case n £ {(0, 0), (1, 0), (0, 1), (1, 1), (2, 0), 
(0,2)}. Note that X^n = 1 at all times, by definition. 
We observe excellent agreement between theory and ex- 
periment. The behaviour of the system can be understood 
as follows. Hydrophilic monomers, with A = 1, prefer lo- 
cal arrangements where contact with the solvent is max- 
imised, i.e. where G {0, 1} (since the maximum number 
of molecules per site is ric — 2). Hydrophobic monomers 
prefer the opposite, i.e. arrangement with n_ e {0,2}. 



9.0 T- 




0.0 10.0 20.0 30.0 40.0 



t 



9.0 p 




0.0 <- ' 1 1 1 1 1 ' 1 

0.0 10.0 20.0 30.0 40.0 



t 

Fig. 2. Evolution of the (co-)variances a± — {n±) and 
p — (n+n_) — (n+)(n_), for the Random Hydrophilic Model 
(RHM), with A' = 2 (A e {-1,1}), Uc = 10, a = i, and 
/3 = 1. Initial conditions: C(2,o) = C(o,2) = 0.5 (perfectly sep- 
arated monomer species). Upper graphs: theory, i.e. numeri- 
cal solution of ( |ll| ) . Lower graphs: numerical simulations, with 
N = 2000. From top to bottom in both pictures: a^{t) (with 
o-i(O) = 1), p{t) (dotted/dashed, with p(0) = 0) and o-+(i) 
(with alio) = 1). 



Hence the observed evolution of the system towards a state 
where only C(o,o)j C(o,2) and C(i.o) are nonzero. 

Note that as soon as AT = 2, with one hydropho- 
bic and one hydrophilic monomer species (as in the ex- 
ample above), the observables in ( [l3|) simply reduce to 
p[n_(.,ri_] — C(„^ „_), and the (co-)variances of 
reduce to those of the cn- In figure (0) we show a first 
example of the evolution of cr^ = (rij.) and p = {n+ri-), 
for the RHM with /? = 1, ric = 10 and a = |. Here 
A+ = 1 and A_ = —1, and the initial conditions were 
C(2.o) = C(o.2) = 0.5 (i.e. perfectly separated monomer 
species). We observe again a perfect agreement between 
theory and (numerical) experiment. The development of 
a large variance cr_ for the hydrophobic monomers is ex- 
plained by the observation that such monomers promote 
arrangements of the type (n+, n_) = (0, ric) or (n+, n_) = 
(★, 0) (i.e. n_ e {0, ric}, without intermediate values). Hy- 
drophilic monomers, on the other hand, promote local ar- 
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bic monomers is to avoid water molecules, and it is ac- 
ceptable for this to be achieved by sharing sites with hy- 
drophilic monomers. Hence, virtually any proposed move 
is accepted which involves a monomer going to a site with 
an adequate number of hydrophobic monomers, regardless 
of the species of the mobile monomer. As a result we ob- 
serve, in addition to the increase of the variance ct^ of the 
hydrophobic monomers, an increase of both the variance 
of the hydrophilic monomers and of the covariance p. 



5 Statics: An Alternative Solution 




t 



Fig. 3. Evolution of the (co-) variances a± = {n±) and 
p — {n+ri-} — {n+){n_), for the Random Hydrophilic Model 
(RHM), with K ^2 {Xe {-10, 1}, i.e. very strong hydropho- 
bicity). He — 10, a — i, and (3 = 1. Initial conditions: 
C(2,o) = C(o,2) ~ 0.5 (perfectly separated monomer species). 
Upper graphs: theory, i.e. numerical solution of (^l|). Lower 
graphs: numerical simulations, with A'^ = 2000. From top to 
bottom in both pictures: o-'-{t) (with cr^(O) = 1), p{t) (dot- 
ted/dashed with p(0) = 0) and o-+(t) (with cr+(0) = 1). 



In the remainder of this paper we will concentrate on 
transitions between low density (swollen) and high den- 
sity (collapsed) states. We distinguish between the overall 
density = N/p (which will remain finite) and the lo- 
cal density as limited by ric, which is the quantity which 
will be allowed to diverge. In [Q the (exact) equilibrium 
solution of site-disordered models was derived in a 

form which was specifically convenient to deal with many- 
valued or continuous disorder distributions; for large ric, 
however, it is very cumbersome (taking the Uc oo limit 
would be ruled out). Since finding collapsed states requires 
ric —> oo, we will first give an alternative solution method, 
which is complementary to that of [Q in that it is conve- 
nient for large Uc, but inappropriate to deal with many- 
valued disorder distributions. In addition it will enable a 
natural and transparent connection with the dynamical 
solution (|Tl|). We write the partition function for the sys- 
tem (IP) as 



.x=l 



(14) 



with the short-hand 



for X > y 
for X < y 



rangements of the type (7i+,n_) = (1,*), i.e. with a single 
optimal value for n+, hence the low variance it+. The low 
observed correlation p indicates that the various objectives 
can be achieved independently. The situation is quite dif- 
ferent for the system illustrated in figure (H), which only 
differs from the previous one in that the hydrophobicity 
of the hydrophobic monomers has been strengthened from 
A = —1 to A = —10. Now the system will find it energet- 
ically favourable to 'absorb' hydrophilic monomers into 
sites where hydrophobic monomers are in contact with 
solvent (that is, in addition to the system aiming for its 
hydrophobic monomers to be in (r;,_|_,n_) = (0, nc) ar- 
rangements, wherever possible) . The energetic gain of the 
hydrophobic monomers, which results upon having their 
hydrophilic cousins replacing solvent molecules, overrules 
the energetic loss of the hydrophilic ones, when the latter 
are taken away from solvent and put into contact with hy- 
drophobic monomers. The primary aim of the hydropho- 



(and its obvious vectorial generalisation 9x.y = Yii ^xi.yi)- 
We introduce the vectors m — (mi, . . . ,mj^) with rrix G 
{0, . . . , ric} to indicate the total number of monomers found 
at each of the R sites. Similarly we define the vector 
fnc = (nc, ■ ■ ■ ,nc). This allows us to write the system's 
free energy as = —(3^^ logX^m (^ma,mZ[m], where the 
'conditioned' partition functions Z[m\ are given by 



Z\m] 



E 



(15) 



We consider only the scaling regime where the limit ric 
oo is taken after the thermodynamic limit N ^ oo. This 
allows us to solve our model in terms of the familiar sub- 
lattices 7a = {«| Ai = A}, of size and their associated 
local occupation numbers n\^j:{r) — X^ig/a ^^-^i- 
ter are, in turn, also grouped into _R-dimensional vectors 
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according to m\{r) — {n\,x{r), . . . ,n\,ii{r)), which al- 
lows us to write the conditioned partition function as 



X (5, 



where 



M! 



is a combinatorial term, counting the number of ways in 
which M objects can be grouped into R subsets of pre- 
scribed sizes {mi, .., mfl}. We note that limAr^oo |^a|/-^ — 
p\ (the a priori probability to draw a monomer of type A 
randomly). In the thermodynamic limit N —f oo the free 
energy per monomer / = lim^v^oo F/N becomes 

A 

m {m^} 



Introducing integral representations for the (5-constraints 
in the terms Af\i^\[mx], i.e. 6nrn = (27r)~i Jo'^du e*"("-'") , 
leads to a full factorisation over sites, and consequently to 
an expression to be evaluated by saddle-point integration 
(provided K remains finite as — s- cx)). This then allows 
us to write the non-constant part of the free energy per 
monomer as: 



/ = -- extr y <^ alog 



E 



=Ea ^^"^"^''Ea.a' 0(^^^')nAnA/ 



A 



(18) 



Variation with respect to y = {y\} gives the K saddle- 
point equations 



a 



1 n,,6 

Ha".! 



'EA"^gEA«-"--i'5EA 



, 0(A,A')TlAn;^/ 

(19) 



with the normalisation factor 

"■-■Ea"^ ^Ea?^^"^-^^Ea.a''^(^'^')"^"a' 



z = 



Da "a! 



(20) 



The relation of this equilibrium solution ( [I8yi9| ) to physi- 
cal observables and to the previously obtained dynamical 
equations ([III ) becomes clear if one adds a suitable in- 
finitesimal generating term to the Hamiltonian (|^) : 



H{r) H{r) 



R 

-y 



n '^nA,x(r),«A 



We repeat the above derivation for the extended Hamilto- 
nian and use the identity lim,j,^o df/dx = (cn (?")), where 
(. . .) denotes the thermal average in the Boltzmann state, 
and where the observables cnix) are defined in (^. This 
reveals that the numbers n = {nx] in (|lq,[l9| ) are identi- 
cal to the local occupation numbers in (|11|), and that in 
equilibrium 



1 



cn 



y\nx 



(21) 



Z Ha "a! 

Hence the saddle-point equations ( [l9| ) simply state that 
(n^) = fjql a. Furthermore, one easily confirms that equa- 
tion (pl|), together with (p^), indeed gives the stationary 
solution of the macroscopic dynamical laws (pT[). 



6 Two Simple Models in the Limit Uc 



We now turn to the simplest case, where there is only one 
type of monomer in the system, such that a local arrange- 
ment at a given site is fully specified by the number of 
monomers present, i.e. n — s- n. Note also that just one in- 
teraction energy </> = (/)(A, A) will be left in our equations. 
Model types will thus only differ in terms of the inter- 
pretation given to the possible values of </>, e.g. for the 
RHM one has (/) = A (^ > for hydrophilic monomers, 
< for hydrophobic monomers), whereas for the RHM 
one has = eA^ (0 > for mutually repelling monomers, 
< for mutually attracting monomers). In this section 
we will for simplicity use the RHM terminology, and re- 
fer to > and </) < as hydrophilic and hydrophobic 
systems, respectively. Equation ( pl| ) now becomes 



m— m! 



(22) 



where Cn gives the fraction of sites with n monomers, and 
where y is to be solved from the non- linear equation (p^, 
which reduces to 

1 ^ 



a 



Tl=0 



n Cn 



This equation can also be written as 
- = |-logZ(y) Z{y) = Y^- 



(23) 



(24) 



I) 



The associated value of the free energy per monomer (|li^ 
is given by 

-/3f = alogZ{y)-y (25) 
We will next study the solution y of (E3j23) and the monomer 



density variance a = 



= iog Z/dy 



the 



different temperature regimes, for large ric- Note that from 
log Z/ dy^ = CT^ > it also follows that the solution of 
(p4|), and thus also of (|2^ , ^ ), is unique. 

In the high-temperature region, both cases (hydrophilic 
and hydrophobic) can be treated on the same footing. Ex- 
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pansion of (p4) in powers of f3 yields 



z{y) = 



(26) 



with Ci — for all other £. For the variance we find 



(33) 



We take the limit 



oo and insert the result into (^ 
which leads to the following quadratic equation for a; = 

1 



The solution gives 

y = -loga 



^/30(i+-; 

2 a 



such that 



1 

a 



1 



aT 



(27) 



(28) 



(29) 



The variance remains finite, and is zero when a~ is an 
integer (in which case a state exists where all sites have 
equal monomer numbers). The system will be in a swollen 
state at any temperature; there are no phase transitions, 
not even for ric — *■ cxd. We now show how this state can 
be extracted from ( |2^ , p3| ) for /9 — > oo. It will be clear 
from ( |2^ ) that a finite a solution requires = 0{T) as 
T ^ 0. Putting y ~ /3<j)x converts ( |23|) into 



1 

a 



E 



ri>0 „! 



Z^n>0 n'. ' 



(34) 



Hence, at finite overall density (i.e. for finite a) and for 
r — !■ oo, the monomer density variance a always remains 
finite (so the system is in a swollen state), due to entropy 
dominance. 

In the low temperature regime, on the other hand, it 
will be clear that the hydrophilic {</) > 0) and the hy- 
drophobic {<j> < 0) models could exhibit qualitatively dif- 
ferent behaviour. More specifically, for models with just 
one monomer type as studied in this section, we can write 
the Hamiltonian (|^) as 



For /3 — !■ oo the two terms n — and n — rix + l will 
dominate the summations in (p^, where G {0, 1,2,.. .} 
such that Ux < X < Ux + ^, and we obtain 



1 n. 
— = lim 

a /3^oo I _(_ 



Upon making the ansatz x = Ua + + log(z)//30, where 
Ua £ {0,1,2,...} such that Ua < a^^ < -I- 1, and 
subsequently taking the /? — > oo limit, we find: 



1 

a 



hmfi 



. z][na + 1] 



H{r) = ^0y^5r,,r, = 7:o:N(f)[^ y^"-x(^) - —] with solution 

2 ' z H ^ — ' a 



^ 1 iV0[a2 + i(i_l)] 
2 a a 



(30) 



lim z — 



1 + hm/3_oo z 



[a ^ - na][l + Ug] 
1 + rin — a^^ 



(35) 



Thus the ground state of the hydrophilic system is the 
one with minimal monomer density variance, and that 
of the hydrophobic system is the one with maximal vari- 
ance. This allows us to determine the ground states with- 
out effort, which provides a welcome test of our general 
equations ( |2^ , |23| ). Extremisation of by variation of the 
{c„}, with two Lagrange parameters for the constraints 
En '''i ~ ^ ^^'^ En ~ Q^^^j shows that an extremum 
will have at most two n with Cn > 0. However, since ac- 
cording to the theory we are not allowed to vary the {cn} 
at will (they must be of the form (p^), any extremum 
thus obtained must be shown to correspond to a specific 
value of the order parameter y in (|2^). 

We first analyse the ground state of the hydrophilic 
model, for ric oo. In order to minimize the energy 
the system must minimise the sum J2i<j^ri,rj, so the 
monomers must spread out maximally over the available 
sites (i.e. the distribution {cn} must be as narrow as possi- 
ble), given the constraints X^n ^n = 1 and J2n '^"^n ~ a^^. 
We define Ua & {0, 1,2,.. .} such that Ua < oT^ < Ua + l. 
The minimum energy arrangement of the {c„} can now 
be written as 



This result is always positive and finite, so our ansatz 
indeed solves the saddle-point equation. We can also cal- 
culate the zero temperature limit of the fraction Cn (p3) 
of sites with n particles. Inserting y = /3(j)[na -I- ^] -f log(z) 
into (p2|) gives 



z e 



ce 



Z_-/rn>0 m! 



> m— o — TT-o 



(36) 



Upon taking the /3 — > oo limit, using (|35[), this reduces to 
the expressions ( pl| , ^2[ ) found earlier by direct variation of 
the {c„}. 

In contrast, the hydrophobic model must maximize the 
variance in order to minimize the energy (see (^0|)), and 
for 4> < the system will consequently compactify (as we 
shall see) at any finite temperature, such that only Cn^ 
and Co are non-zero: 



1 = 0: Co l-(anc)"^ 



(37) 
(38) 



ce 



(31) 
(32) 



with C£ = for all other £. The variance then equals, in 
leading order as —» oo. 



Tic/ a — Xj o? 



(39) 



8 



H. Chakravorty, J. van Mourik, A.C.C. Coolen: Solvable Lattice Gas Models of Random Heteropolymers II 



and consequently diverges. Apparently, for Uc oo the 
equilibrium state of the hydrophobic system is one with 
just two types of site: those with Uc particles and those 
with no particles (i.e. it is a maximally compact state), 
for any temperature. We will now study the ric oo be- 
haviour of the saddle-point equation (p3|), for </> < and 
arbitrary (3. We make the ansatz 



1 



(3(j)nc 



1 



log(a7^) -t- \og{nc) + X (40) 



where lim< 



Uc — '■OO ■ 



0. This converts (p4) into a ^ 



9 log Z{x)/dx, with, upon application of Stirling's formula 
to the term n — ric in Z{y), and modulo contributions to 
the exponents which vanish as ric ^ oo: 



Z{x) = 1 -f 

ric-l 

+ E 

n=l 



aricvi-K 



!{...) 



1 r 



^2;-l+i/3|0|(n-nc)--^ log(ayri7)+log(nc) 



One distinguishes two qualitatively different types of terms 
in the above summation. Firstly, those with n <^ ric oi 
with n = Oiric) (provided n <C ric) give a vanishing cu- 
mulative contribution to Z{x) due to the appearance of 
\(3\(j)\{n — ric) in the exponent. Secondly, those terms with 
n = ric — k, with finite fc > 0, each give according to Stir- 
ling's formula and modulo contributions to the exponents 
which vanish as ^ oo: 



1 



^2;-l + i/3|0|(n-nc)-;7- log(Qyn7)+log(nc) 



aricV^TT 



■^naP\<t>\k-k(x-l-^f3\4i\k)+0{...) 



Hence also the second type of terms in Z{x) are exponen- 
tially damped as ric ~* oo, relative to the two surviving 
terms n = and n — ric, and hence 



Z{x) = 1 + 



[l + 0{. 



aUcV^Ti 

Thus for ric ^ oo the saddle-point equation reduces to 



(41) 



1 



lim 



ricC-^ 



and hence x = n^^ log ^/2^T -\- . . . as ric ^ oo, confirming 
the correctness o f our a nsatz. This gives for the fractions 

maximisation ^ 



(22) the state (|3^ , |38| ) which resulted earlier from direct 
,ximisation of by variation of {c„}. 



7 The Transitions to Compact States 

We have observed in the simple single-species models of 
the previous section that in the hydrophilic case there are 
no phase transitions for ric ^ oo (since here the entropic 
and energetic forces have the same objective: a swollen 



state), whereas in the hydrophobic case the system will 
for ric —> oo be found in a maximally compact state at any 
temperature (i.e. there is a phase transition at T = oo). 
Since the latter behaviour is strongly linked to the ric oo 
limit, one should expect that proper phase transitions in 
the hydrophobic system can be identified in terms of a 
suitably re-scaled nc-dependent temperature (or, equiva- 
lently, a suitably re-scaled interaction energy 0). It turns 
out the the proper scaling is obtained upon transforming 
according to 

T-f 101^2/ log(nj), (42) 

(where we have absorbed the interaction energy into the 
re-scaled temperature), together with /3 = 1/T. For ric — > 
oo we observed that the fraction of sites with ric parti- 
cles vanished, even in the maximally compact state, where 
Cric = (c^c)"^)- Hence we switch in this section from in- 
specting the fraction c„ ( ^ ) of sites with a given number n 
of particles, to an alternative representation of the macro- 
scopic system state in terms of the fraction of particles 
/„ = anc„ which find themselves at a site filled with n 
particles (in the maximally compact state one has = 
for all n < ric and = 1). Note that /o = 0, by defini- 
tion. We introduce the short-hand 



h{n,y) 



jdn"^ log(ne!) 
2^2 



log(n!) + ny (43) 



The saddle-point equation (|23|), from which to solve y, can 
now be written as 



^ ^ fn — Ij fn — Y~>" 



:,h{n,y) 



(44) 



Equation (^4|) looks deceivingly simple. However, it turns 
out that one cannot simply transform the previous so- 
lution y (^) of the [3 — C(?^c) regime by writing [3 in 
terms of the re-scaled (3, due to the non-commutin g of the 
various operations involved. Numerical solution of ( [44| ) in- 
dicates that for ric ^ oo we can neglect all those terms n 
in the summations with \ <^ n < ric, relative to n — ric, 
and that terms with n ~ C(l) are generally of the same 
order of magnitude as that of n = 1. One can show that 



this property can be derived, from (44), as an immediate 
consequence of assuming that the leading orders of /i and 
fn^ scale in the same way with ric as ric ~^ oo. The latter, 
in turn, is equivalent to saying y = (1 — kP) log »-c + u 
with u finite for nc — > oo, exactly as in (Eo) (after car- 
rying out the transformation from (3 to 0). As a result, 
and in contrast to the scaling regime (3 = O(n^), even 
for nc ^ oo equation ( ^ ) remains transcendental, and it 
generally resists analytical solution. 

A systematic numerical analysis of (Q) reveals the fol- 
lowing picture. At high temperatures, /? = 0, the system 
is in a swollen phase, which is defined by the condition 
lim„^_,oo /ric = (the fraction of monomers found at a 
maximally filled site is zero) . Upon increasing (3 one then 
finds a phase transition at /3f(a) to a partially compact 
state, where < lim,i^^oo fn^ < 1. A further increase of 
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6.0 p 

5.0- compact = 1) 



4.0 -^^^ 

3 partially compact (0</„^<l) 



2.0 

- swollen (/„,=0) 



0.0 1.0 2.0 3.0 4.0 5.0 

a 

Fig. 4. Phase diagram of the hydrophobic model (with a sin- 
gle monomer species), as obtained by solving the saddle-point 
equation numerically (with = 2.10^). The re-scaled 

inverse temperature /3 is defined in (^). 

P leads to a growth of the fraction of monomers in maxi- 
mally filled sites, until at a second transition temperature 
/?|(a) the system enters the fully compact phase, where 
lim„^_»oo /ric — 1- For the special case /3 = 2 it turns out 
that ( ^ ) can be solved analytically: for a > 1 one finds in 
leading order y = — loga and lim„^_>oo fn^ = 0, for a < 1 
one has y = 0{n~^ logUc) and < lim„^.^oo /n, < 1, 
which confirms the picture obtained by numerical solu- 
tion (see below). 

The resulting phase diagram, obtained by extracting 
the two transition temperatures /3J 2 i'^) from a numerical 
solution of (|4^), is drawn in figure (||), and is confirmed 
qualitatively by numerical simulations (see the next sec- 
tion). Note that even for the relatively high value Uc = 
2 X 10^ as used in the numerical evaluation, the line P^^i'^) 
can only be qualitatively correct, since in order to probe 
truly the ric c>o regime we need log(nc) — > 00, according 
to the saddle-point equation. Note, furthermore, that in 
the limit a —^ the system must be fully compact, so that 
both f3i and /Sf must tend to 0. Finally, for a 1 our nu- 
merical solution shows that both transition temperatures 
grow as Pi,2{a) ~ log(a). 



8 Numerical Simulations for High Uc 

We have carried out numerical simulations to verify the 
theoretical equilibrium results of the last section, and the 
kinetics of the approach to equilibrium in the ric — > 00 
regime; since phase transitions only occur for ^ < (the 
hydrophobic system) we have restricted our simulations 
accordingly. Due to often excessive equilibration times ob- 
served for large ric (see below), it has been necessary to 
limit our experiments to ric = 100, N — 1000. Hence, 
strong finite size effects could be anticipated, since the 



1.0 



0.2 



- 


J 


1 




/ 

/ 








• • Collapsed phase 

Swollen phase 









) I ——i i e ' o 00 1 !. ' i ' i ' 1 

0.0 2.0 4.0 6.0 8.0 10.0 



/3 

Fig. 5. Fraction of monomers in a maximally filled site, 
measured in numerical simulations of the single-species hy- 
drophobic system ((^ < 0), as a function of the re-scaled in- 
verse temperature. Simulation parameters: N = R — 1000 
(so a = 1), ric = 100, and relaxation time t — 10® iter- 
ations/monomer. Initial conditions: either the compact state 
(solid), or the swollen state (dashed). 

theory requires N —>■ 00, ric —>■ 00 and ric/N 0. Nev- 
ertheless, the numerical simulations were found to agree 
remarkably well with the theory. 

Figure ^ shows the fraction of monomers in a max- 
imally filled site, as a function of the re-scaled inverse 
temperature /3, after 10® updates per monomer. The solid 
curve shows the observed value following initialisation of 
the system in a maximally compact (collapsed) state; the 
dashed curve shows the observed following initialisa- 
tion in the fully swollen state. The hysterisis-type graph 
which results from these measurements (the dependence 
on intialisation) refiects extremely slow equilibration of 
the system towards more compact states, due to entropic 
barriers (similar to those in e.g. |lj]). The full curve (which 
describes a process of de-compactification) does not suffer 
from this problem, and therefore gives the best experimen- 
tal estimate for the true locations of the transitions. This 
latter curve predicts (for a = 1.0) that — 1.7 and 

/32(1) — 4.1, to be compared to the theoretical predictions 
/3f(l) ~ 2.0 and /3f ~ 4.2 (see figure |). Given the antici- 
pated finite size effects, these values are remarkably close. 
The non-monotonicity of some of the experimental curves 
are surely finite ric effects; for finite ric a maximally filled 
site can indeed temporarily loose one or more monomers, 
whereas the probability for this to happen vanishes for 
ric oo- 

We now turn to the entropy barriers which cause the 
slow equilibration towards compact states. At sufficiently 
low temperatures, starting from a swollen (i.e. high tem- 
perature) state, monomers initially hop around until at 
some site a critical (temperature dependent) monomer 
number is obtained. After this has happened, this site fills 
up very rapidly, since its monomers can no longer escape 
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Fig. 6. Variance — (n^) — (n)^ (proportional to the energy) 
as a function of time, measured in numerical simulations of 
a hydrophobic system as in the previous figure, at /3 = 3.8, 
following fully swollen initial conditions. 



(for energetic reasons). All monomers in fully filled sites 
are thus more or less immobile, and the remaining (mo- 
bile) monomers find themselves at a reduced effective den- 
sity Pcff = (1 — fn^)/ict — Cn^)- Hence, as increases, for 
the mobile monomers it will become increasingly harder 
to 'find' other monomers to form a new site with a critical 
monomer number. This entropic effect is illustrated by the 
simulation shown in in figure |[ It is characterised by sud- 
den jumps in the variance, corresponding to a site reaching 
the critical monomer number and filling up rapidly, with 
increasingly long plateaus in between, where the remain- 
der of mobile monomers hop around at a reduced effective 
density. The actual process of the rapid filling of those 
sites which a some stage reach the critical monomer num- 
ber, is illustrated in figure]^. It shows, as snapshots of nu- 
merical simulations at three intermediate (relatively close) 
times, the distribution {/„} for a hydrophobic system 
as in the previous experiments which has been quenched 
from a maximally swollen state to a lower temperature, 
well within the fully compact phase. Since the equilib- 
rium state in that phase corresponds to — 1, even at 
the largest observation time shown, = 10'', the system 
clearly has not yet equilibrated. 

Another interesting aspect of the equilibration dynam- 
ics is the fact that compactification is found to take longer 
for a system which is not re-initialised at each new relax- 
ation temperature, compared to a system that is. This 
is illustrated in figure @. Apparently, in a quench (re- 
initialisation to a swollen state), several sites reach the 
critical monomer density simultaneously, while the effec- 
tive density is still reasonably high. Note that at short 
times one can indeed observe several peaks in figure ^. 
On the other hand, in a system which is not re-initialised, 
the effective monomer density is already reduced when 
the temperature is changed and a new relaxation com- 
mences. As a result, the entropic barrier to be overcome 
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Fig. 7. The occupation fractions /„ as observed in a single 
numerical simulation at times ti = 10^, t2 —2 10^, ta — 10'*, 
for a hydrophobic system with N — R = 1000, ric = 100, and 
/3 — 8.0. It illustrates the formation, within a relatively small 
number of iteration steps, of a maximally filled site. 

before a critical monomer density at any site is achieved, is 
larger in the system which is not re-initialised than in the 
re-initialised system, and hysteresis effects are therefore 
stronger. Non-monotonicity is, again, a finite rig effect. 
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Fig. 8. Occupation fractions versus the re-scaled inverse 
temperature /3, as measured at t = 1 10^ in simulations of a 
hydrophobic system with N = R = 1000, ric = 100. Upper 
graph: the system is not re-initialized when starting a new 
simulation at increased (3. Lower graph: for each /3 the system 
is re-initialized to a maximally compact or swollen state. 

9 Discussion 

In this paper we have studied the dynamics and the com- 
pactification phase transitions of a class of simple lattice 
gas models for a random heteropolymer in infinite dimen- 
sions, as introduced in Within this class we have re- 
stricted ourselves to models with site disorder, in particu- 
lar to the Random Hydrophobic Model (RHM). For finite 
ric (the maximum number allowed at a given site), and 
a finite number of monomer species, we have been able 
to solve the dynamics exactly, by deriving exact closed 
macroscopic evolution equations for a suitable set of dy- 
namic order parameters. The solution of the (non-linear) 
macroscopic equations, describing the approach to equilib- 
rium, is shown to be in excellent agreement with numerical 
simulations. 



In order to study the physics of the model also for 
ric oo (a limit which could not easily be taken in the 
equilibrium solution proposed in |^ ) we have presented an 
alternative (but equivalent) equilibrium solution. In the 
latter approach the limit ric oo poses no problems (al- 
though the method would not apply to models with bond 
disorder, in contrast to that of ||l[). Working out our equa- 
tions for models which just a single monomer type shows 
that for Tic — > oo a hydrophilic model will not undergo 
phase transitions at any temperature (as expected), but 
that the hydrophobic system will be in a collapsed state, 
at any finite temperature. In the latter case the transi- 
tion(s) apparently occur at T = oo, as a consequence of 
having taken the ric ^ oo limit. In terms of a suitably re- 
scaled inverse temperature /3, however, two transitions are 
found to occur, at (3f (from a swollen to a partially com- 
pact state) and at (from a partially compact to a fully 
compact state) . Both the existence and location of these 
two transitions are confirmed by numerical simulations, 
within the experimental limits imposed by having finite 
TV and finite ric (the theory requires N oo, ric oo 
with ric/N 0). Especially when following a swollen (i.e. 
high temperature) initial state, equilibration is observed 
to take place via an extremely slow plateau-like dynam- 
ics, which is caused by entropic effects, particularly at the 
critical points and in the partially compact phase. 

A natural extension of the work reported in this paper 
would be to solve the dynamics of the model for the case 
of bond disorder. In statics such models have been solved 
using rephca theory, see ||l|,||,|4[||,P,||,||], and it is to be 
anticipated that the method of path integrals ||l^ would 
have to be used for solving the dynamics, as in e.g. 0. 

Infinite dimensional models are obviously less realistic 
than finite dimensional ones in terms of the spatial prop- 
erties of polymer folding processes. However, our present 
study shows that, due to the absence of any need for ap- 
proximations (in contrast to their finite dimensional coun- 
terparts), they do allow for a fairly complete and fully 
exact quantitative understanding of the energetic and en- 
tropic interplay between hydrophobic and hydrophilic ef- 
fects, including both statics and dynamics. 
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